DC and AC Josephson Effect in a Superconductor-Luttinger Liquid-Superconductor 
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We calculate both the DC and the AC Josephson current through a one- dimensional system of 
interacting electrons, connected to two superconductors by tunnel junctions. We treat the (repulsive) 
Coulomb interaction in the framework of the one-channel, spin-1/2 Luttinger model. 
The Josephson current is obtained for two geometries of experimental relevance: a quantum wire 
and a ring. At zero temperature, the critical current is found to decay algebraically with increasing 
distance d between the junctions. The decay is characterized by an exponent which depends on the 
strength of the interaction. At finite temperatures T, lower than the superconducting transition 
temperature Tc, there is a crossover from algebraic to exponential decay of the critical current as a 
function of d, at a distance of the order of hvp/ksT. Moreover, the dependence of critical current 
on temperature shows non-monotonic behavior. 

If the Luttinger liquid is confined to a ring of circumference L, coupled capacitively to a gate voltage 
and threaded by a magnetic flux, the Josephson current shows remarkable parity effects under the 
variation of these parameters. For some values of the gate voltage and applied flux, the ring acts 
as a TT-junction. These features are robust against thermal fluctuations up to temperatures on the 
order of hvp/kBL. 

For the wire-geometry, we have also studied the AC- Josephson effect. The amplitude and the phase 
of the time- dependent Josephson current are affected by electron-electron interactions. Specifically, 
the amplitude shows pronounced oscillations as a function of the bias voltage due to the difference 
between the velocities of spin and charge excitations in the Luttinger liquid. Therefore, the AC 
Josephson effect can be used as a tool for the observation of spin-charge separation. 



PACS numbers: 74.50 -fr, 72.15 Nj 
I. INTRODUCTION 

Due to the recent development of superconductor- 
semiconductor (S-Sc) integration technology it has be- 
come possible to observe the transport of Cooper pairs 
through S-Sc mesoscopic interfaces. Examples are the 
supercurrent through a two-dimensional electron gas 
(2DEG) with Nb contacts (S-Sc-S junction) 0] and ex- 
cess low-voltage conductance due to Andreev scattering 
in Nb-InGaAs (S-Sc) junctions Hj. The transfer of single 
electrons through the interface between a semiconductor 
and a superconductor with energy gap A is exponen- 
tially suppressed at low temperatures and bias voltages 
kBT,eV ^ A (e is the electron charge). Instead, elec- 
trons will be transferred in pairs through the interface, 
a phenomenon known as Andreev-reflection It has 
been realized only recently, that the phase coherence be- 
tween the two electrons involved in this process could give 
rise to distinct signatures in the transport properties of 
mesoscopic S-Sc-S and S-Sc systems Q,^. 

If the normal (Sc) region is free of disorder, the prop- 
agation of electrons is ballistic. Phase coherence be- 
tween the two electrons is maintained over the length 
Lcor = hvp/'raax {ksT, eV}, where vp is the Fermi ve- 
locity. In this regime, the critical current Ic through a 



short and narrow constriction in a high-mobility non- 
interacting 2DEG, connected to two superconductors 
should be quantized each propagating mode con- 
tributes an amount eA/h to the critical current. 

In the presence of disorder in the normal region, the 
motion of the two electrons will be diffusive. Like in 
disordered metals, the phase coherence between the two 
electrons is limited by the correlation length Lcor = 
yJfiD/m.&yi {ksT, eV}, where D is the diffusion constant. 
For instance, the excess low-voltage conductance in S-Sc 
junctions can be explained in terms of constructive 
interference occurring over this length scale between the 
two electrons incident on the S-Sc interface [@j8|. 

In these examples, electron-electron interactions are 
neglected. It is well known, though, that they may have 
a strong influence on the transport properties of meso- 
scopic systems. In general the interactions modify the 
phase-coherence length Lcor, which poses limitations on 
the abovementioned mesoscopic effects. In specific cases 
the effects of electron-electron interactions will strongly 
depend on the layout of the system under consideration. 

For example, the interactions will modify the critical 
current Ic through a normal metallic slab sandwiched be- 
tween two superconductors Q . If the coupling between 
normal metal and superconductors is weak (tunneling 
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regime) and the size of the slab (and hence its electric 
capacitance C) is small, a phenomenological capacitive 
model [0 can be used to describe the effect of interac- 
tions. As a result the critical current shows strong reso- 
nant dependence on the electro-chemical potential of the 
slab dependence has different character for Ec < A and 
Ec > A, Ec = e^/2C being the charging energy. On the 
other hand, if the normal metal and the superconductor 
are well coupled (regime of Andreev reflection), electron- 
electron interactions will modify the results obtained in 
Ref. 1^ in quite a different fashion. A perturbative treat- 
ment of the interactions |jl^ shows that an additional 
supercurrent through the slab arises, whose sign depends 
on the nature of the interactions in the slab (attractive 
or repulsive), and whose phase-dependence has period tt 
(rather than 2tt in the non- interacting case). 

If instead of a metal a low-dimensional Sc nanos- 
tructurc with a small electron concentration is consid- 
ered, the abovementioned descriptions of the electron- 
electron interactions are no longer sufficient. In one- 
dimensional (ID) systems the Coulomb interactions can 
not be treated as a weak perturbation. As a result a 
non-perturbative, microscopic treatment of interactions 
is required. For ID systems this can be done in the frame- 
work of the Luttinger model Interactions have a 
drastic consequence: there are no fermionic quasiparti- 
cle excitations. Instead, the low energy excitations of 
the system consist of independent long- wavelength oscil- 
lations of the charge and spin density, which propagate 
with different velocities. The density of states has power- 
law asymptotics at low energies and the transport prop- 
erties cannot be described in terms of the conventional 
Fermi-liquid approach. For a quantum wire with an ar- 
bitrarily small barrier this leads to a complete supression 
of transport at low energies . 

Another interesting feature arises in ID interacting 
systems of a finite size. For a Luttinger liquid confined 
to a ring. Loss found remarkable parity effects |l^] 
for the persistent currents. He used the concept of Hal- 
dane's topological excitations |jT^, extending the pre- 
vious work of Byers and Yang for non-interacting elec- 
trons in a ring |19| . Depending on the parity of the to- 
tal number of electrons on the ring, the ground state is 
either diamagnetic or paramagnetic. For spin-1/2 elec- 
trons an additional sensitivity on the electron number 
modulo 4 has been foimd ||20|| . Experimental evidence 
for Luttinger liquid behavior in Sc nanostructures has 
been found recently. The dispersion of separate spin 
and charge excitations in GaAs/AlGaAs quantum wires 
has been measured with resonant inelastic light scatter- 
ing pl| ] . Transport measurements on quantum wires have 
revealed power-law dependence of the conductance as a 
function of temperature p^ . 

In view of this we expect that electron-electron inter- 
actions may well have drastic, observable consequences 
in systems which consist of low-dimensional Sc nanos- 



tructures connected to superconductors. In this paper, 
we will study the Josephson current through a Lut- 
tinger liquid [ p3|j2^ ]. Specifically, we consider two ge- 
ometries which can be realized experimentally: a long 
wire with contacts to two superconductors at a distance 
d (see Fig. ||a) or a ring shaped Luttinger liquid shown 
in Fig. |l|b. In both cases the ID electron liquid is con- 
nected to the superconducting electrodes by tunnel junc- 
tions. This is an interesting system from various points 
of view. First of all, it enables one to study in a mi- 
croscopic way how the Coulomb interaction influences 
the phase-coherent propagation of two electrons through 
a ID normal region ||2^. Secondly, various aspects of 
transport in mesoscopic systems (parity effects and in- 
terference combined with electron-electron interactions) 
and their interplay can be enlightened in such a device. 
Finally, since the Josephson effect is a ground state prop- 
erty, the Josephson current can be used as a tool to probe 
the ground state of an interacting electron system. In 
particular, for the ring geometry in the presence of an 
Aharonov-Bohm flux, the various possible ground-state 
configurations can be determined by studying flux and 
gate voltage dependence of the critical current. 

The paper is organized as follows. In Section 2 we 
briefly review the properties of the spin-1/2 Luttinger 
model. In Section 3 the general formalism for the DC- 
Josephson effect is presented. The DC- Josephson cur- 
rent is obtained by evaluating the contribution to the 
free energy which depends on the difference of the su- 
perconducting phases. Starting from the general expres- 
sion (|l2|) we then consider various interesting limiting 
situations. In Section 4, the wire geometry is considered 
(see Fig. la). The critical current decays as a power 
of the distance d between the contacts. The exponent 
depends on the interaction strength. We distinguish two 
cases in which the characteristic energy hvp/d for the ID 
system is either much smaller or much larger than the su- 
perconducting gap A. For the ring geometry (Section 5, 
see Fig. lb), we focus on the dependence of the critical 
current on the applied gate voltage and/or flux. Both 
Section 4 and 5 contain a discussion of the effect of finite 
temperatures. Section 6 is devoted to the AC-Josephson 
effect. In this case the imaginary time approach of the 
previous sections is inadequate and we will use a real 
time formulation. The amplitude of the AC component 
is found to show oscillations as a function of voltage due 
to spin-charge separation. In the last section we present 
the conclusions. 



II. THE SPIN-1/2 LUTTINGER LIQUID 

We start the description of the model we use by re- 
viewing the theory of ID interacting spin-1/2 fermions 
(throughout, we use ?i = fcs = 1)- The long- wavelength 
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behaviour of such a quantum wire of length L is governed 
by the Hamiltonian 



Hr 



L/2 

/dx \ ^ 



-L/2 



9i 



It is written as a sum of the contributions from the 
spin (j — a) and charge (j = p) degrees of freedom. 
The parameters gj denote the interaction strengths and 
Vj = 2vp/gj the velocities of spin and charge excita- 
tions iQ. The parameters gj can be determined once 
one defines an appriopriate microscopic Hamiltonian (e.g. 
the Hubbard model); an approximate form for the spin- 
less case has been given in Refs. In this paper 
we will neglect backscattering (via Umklapp or impu- 
rity scattering) and restrict ourselves to repulsive, spin- 
independent interactions. As a result we have g^- = 2 and 

V„ = VF 0. 

We also introduced bosonic fields 0j and dj. They are 
related to the fields 4)s — (t>p + s(f>a and 6s = 9p + sO^ for 
spin up (s — +1) and down (s = —1) fermions. These 
fields obey the commutation relation [0s(x), (x')] = 
(i7r/2) sign(x' — x)5s,s' ■ The fermionic field operators ^ 
can be expressed in terms of the spin and charge degrees 
of freedom (l|Jl|: 

n,odd 

X exp{m[6'p + s6'cr]} exp {i[0p + scj)a]} , (2) 

where kp is the Fermi- wave vector and po = Nq/ L is the 
average electron density per spin direction. The num- 
ber Nq determines the linearization point of the original 
electron spectrum, hp = ttNq/L |16|. 

If the wire is closed to form a loop, the periodic con- 
dition 



The boundary condition (||) gives rise to the topologi- 
cal excitations Mj and Jj for the spin and charge degrees 
of freedom. They are related to the usual topological ex- 
citations for fermions with spin s: Ms = {l/2)[Mp+sMa-] 
(1) and Js = (l/2)[Jp -I- sJ„]. Physically, the number Ms 
denotes the number of excess electrons in the Luttinger 
liquid (LL) with spin s in addition to Nq. The number 
Js is the number of current quanta evp/L, carried by 
electrons with spin s. Here, vp = ■RpQ/m is the Fermi 
velocity, with m being the electron mass. A net current 
JpCVp/L flows through the quantum wire if there is an 
imbalance between the number electrons moving to the 
right and to the left. Using the boundary condition (|3|) 
one obtains topological constraints for Ms and Js ||l6|,|l8[ , 
which lead to the following constraints for Mj and Jj-. 
(i) The topological numbers Mj and Jj are either simul- 
taneously even or simultaneously odd; 



i>LA^ + L) = (-i)^"'j'i,,(x). 



(3) 



should be imposed on the Fermi operators (|^) . The fields 
9 and (/) can then be decomposed in terms of bosonic fields 
9 and (j) and topological excitations M,18[ : 



,{x) = 9j{x) + 9"^ + ttMj{x/2L), 



(4a) 



(x) = 4>j (x) + <j)"j + TT Jj [{x + L/2)/2L] . (4b) 

Here, 0j and are given by 
1/2 




3ign(g)e'«-(6]_^ + 6,-_,), (5a) 



(5b) 



where bj^q, &j ^ are Bose operators. 



(ii) when A^o is odd the sum Mp±Mcp 



'Jp±Jc 



takes values 



—4, 0, 4, when Nq is even the sum Mp±Mrj + Jp±Jcr 
takes values ...,—6,-2,2,6,.... 

An Aharonov-Bohm flux $ threading the loop couples 
to the net current, characterized by the topological num- 
ber Jp of the field (f)p. The flux can be incorporated into 
the Hamiltonian (|l|) by changing 

V(/., ^ - (27r/L)<5,- p/$, (6) 

where /$ = ^/^o is the flux frustration and $0 is the 
flux quantum h/e. 

Since the LL is brought into the contact with par- 
ticle reservoirs (superconductors) kept at fixed electro- 
chemical potential /i, the number of particles (character- 
ized by the topological number Mp of the field 9p) should 
be coupled to p. This can be achieved by replacing 



V9, - i2n/L)6,^pf, 



(7) 



in the Hamiltonian (0). The parameter fp = 
{gpL/4:TrVp)Afi is related to the the difference Ap be- 
tween /i and the Fermi-encrgy Ep — kp /2m of the quan- 
tum wire, corresponding to the linearization point. Gen- 
erally, the reference point A/i = is defined from the 
requirement that for $ = there are 2A^o electrons in 
the ground state and the energies to add/remove elec- 
trons to/from the system are equal. The difference A/i 
can be controlled, e.g., by a gate voltage. 

Using Eqs. (§), (11), (^), and 1^ one concludes 
that the Hamiltonian can be decomposed into non-zero 
modes and topological excitations: 



Trvj_ 
4L 



%{J, - A5,.pU f + -{M, - 4S,,pU f 
^ 9j 



(8) 



Since this Hamiltonian is quadratic in the Bose opera- 
tors, it is possible to obtain all the correlation functions 
exactly. 
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III. DC-JOSEPHSON EFFECT 

Both systems depicted in Fig. || can be described by 
the Hamiltonian 



Fs^{0■,T-T') = (TAs^-{d,T)^S^,+ {d,T') 



Si 



H — Hsi + Hs2 + Hl + Ht = -ffo + Ht- 



(9) 



Here, Hsi, Hs2 are the BCS-Hamiltonians for the bulk 
superconductors kept at constant phase difference x — 
Xsi ~ XS2- For simpHcity we assume equal magnitude A 
for both order parameters. The tunneling between the 
superconductors and the ID electron system is described 
by It is assumed to occur through two tunnel junc- 
tions at the points a; = and x = d, 

i?T = ^Ti*^i,(x = 0)^'L,.(x = 0) 

S 

+T2^S2A^ = d)^LAx = d) + (h.c). (10) 

The constant tunnel matrix elements Ti_2 can be related 
to the tunnel conductances G1.2 of the junctions, Gi = 
ATre^NLiO)Ni{0)\Ti\^, where Nl{0) = I/ttwf, and N, is 
the normal density of states in the superconductors {i = 
1,2). 

The stationary Josephson effect can be obtained by 
evaluating the phase-dependent part of the free energy 
J^ix)- The Josephson current is then given by 



-2edT/dx. 



(11) 



We expand T = -(l//3)lnZ, where Z = Tr exp{-PH} 
and /? ~ 1/r, in powers of the tunneling Hamiltonian 
Ht, using standard imaginary-time perturbation the- 
ory . The lowest order phase-dependent contribution 
arises in 4"* order. Using Eq. (|l^) we see that there are 
24 contributions to the phase-dependent part of T: 

Ti T2 TA 

T{x) = ~^ j dri j dT2 J drs J dr^ 
0000 
X {4i(0; n - T2)T^Tl[^\0, d- n, Ti) 

y-{nfFs2{Q\T^-Ti) 
+4^(0; n - Ti)Tln^l^\Q, d- n, T4) 

x(T2*)2i^S2(0;T2-T3) + (h.C.)} 

+22 similar terms. (12) 

This result has a clear physical meaning, see Fig. |. The 
Josephson effect consists of processes in which a Cooper 
pair tunnels from superconductor S2 into the LL with an 
amplitude (12*)^. After propagation through the LL, it 
tunnels into superconductor SI with an amplitude (Ti)^. 
The hermitian conjugate terms describe processes in the 
opposite direction. The propagation in the superconduc- 
tors is described by the anomalous Green's function 



7riV(0) 
~1~ 



(13) 



where {■■■)si indicates an average with respect to Hsi- 
Propagation through the LL is determined by the 
Cooperon propagator Hi(0, d; ti, T4). These 24 terms 
are obtained by considering all possible time-ordered 
pairs of tunneling events , Tj (r^ < Tj) at x = and 
X — d together with all possible spin configurations. 
However, which of these terms are important depends 
on the relation between the characteristic energy vp/d 
for the ID system and the superconducting gap A [p8| . 

If the distance between the contacts is large, vp/d <C 
A, a generic process consists of fast tunneling of two elec- 
trons from the superconductor into the ID system and 
their slow propagation through the LL. Such a process is 
illustrated by the first term in Eq. (p^. Here, 

n(:')(o,d;n,...,r4) 

= (#L,+ (0,Tl)*i^_(0,T2)^'[_+(d,T3)*l__(d,T4)), (14) 

where |ti — ra | ^ d/vp > In — ^2 1 ~ |t3 — r4 1 ^ 1/ A (see 
Fig. |a). The average is taken over equilibrium fluctua- 
tions in the LL (described by Hl)- The other relevant 
processes come from diagrams which are obtained from 
the one in Fig. ^ja by means of particle-hole conjugation 
and by changing the time ordering. 

In the opposite limit vp/d 3> A, diagrams of the type 
depicted in Fig. |^a are no longer relevant. Instead, one 
should consider fast and independent propagation of two 
electrons through the LL and slow tunneling between S 
and LL. This is illustrated by the second term in Eq. (|lT 
(see also Fig. |^), where 



ni''(0,d;ri. 



, T4) 



= {^^L.+iO,ri)^U{d,T2)^L,-iO,Ts)^l_id,n)j, (15) 

with In - T2I |t3 - r4| ~ d/vp < |ti - T3I ^ 1/A. 
Also in this case the other relevant processes can be ob- 
tained from the one in Fig. |b by means of particle-hole 
conjugation and by changing the time ordering. 

The direct evaluation of averages like ([l|), (|l|) with 
the help of bosonized field operators like (g) is tedious 
but straightforward. The resulting expressions can be 
simplified further in the two limiting cases vp/d <C A 
and vp/d^ A, which contain all the important physics 
of the problem. 



IV. DC-JOSEPHSON CURRENT THROUGH A 
QUANTUM WIRE 

We first turn to the geometry depicted in Fig. |l|a. It 
consists of a quantum wire of length L ^ 00 connected 
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to two superconductors by tunnel junctions separated by 
a distance d. The topological excitations play no role in 
this case (their energy is vanishingly small) and the wire 
is described by the non-zero modes only (first term in 
Eq.(|)). 

a. The case vp/d <^ A. 

The expression for the phase-dependent part of the 
free energy contains four terms of the type of the first 
term in Eq. ( p^ (see Fig. |2^). In this case at low tem- 
peratures T ^ A one can approximate the anomalous 
Green functions (^3|) by (5-functions in time. This fixes 
equal time arguments ti — T2, = T4 in n^"\ The re- 
maining integration should be performed over the time 
T — Ti — T3. The dominant contribution to the the inte- 
gral comes from the terms with n = ±1 in Eq. (||). As a 
result, the Josephson current ( [III ) through the quantum 
wire is given by 



surx, 



with a temperature-dependent critical current 



Airevp G1G2 



(4e2)2 



(16) 



(17) 



where 



kpd 



2/3P-1 



dx y-r 

2^ J^J^ 



2^ 



cosh{2TTd/vjP) — cos {2tt dx / (3vf) 



(18) 



(with — 2 and — vp)- 

In the non-interacting case, at zero temperature, 
Fw\o) = 1- The Josephson current decreases as 1/d 
with increasing distance between the tunnel junctions. 
This is related to the fact that the density of Cooper pairs 
in the LL decays in space away from each junction. Hence 
the overlap of the macroscopic wave functions of the two 
superconductors, which is responsible for the Josephson 
effect, is suppressed. Repulsive interactions in the wire 
make the Josephson effect vanish more rapidly with the 
distance between the superconductors, 

li'^^O) Oil/ {kpdf . 

The electron liquid acquires an additional stiffness 
against density fluctuations, hence the tunneling between 
S and LL is suppressed. This fact provides an a posteri- 
ori justification of our use of perturbation theory when 
treating electron tunneling in the presence of repulsive 
interactions. 

We consider now the temperature dependence of the 
critical current. For non-interacting electrons, gp = 2, 
the critical current can be calculated explicitly: 



/i"^(r) _ 2nTd 



cosh^{2TTTd/vp) - 1 



(19) 



At low temperatures, T <C vp/2Trd, the critical current is 
suppressed below its zero temperature value in a power 
law fashion 



^ 2 /irTd 
3 \ vp 



(20) 



In the high temperature regime, T <C vp/2-Kd^ the decay 
is exponential 



ir{T) ^ VSirTd 



Vp 



exp(—2TrTd/vp) 



It is possible to obtain analytical results also in the 
interacting case. In particular, for weak interaction 2 — 
gp <^ 2 and low temperatures T <^ vp/2'Kd the critical 
current behaves as 



I^"\t) ^ 2-gp /2TTTd 
&\0) 4. \ VP 



2 fnTdy 



Vp 



, (21) 



where we dropped terms 0[{2- gp)"^ + {2- gp){Td/vpf]. 
This result can be interpreted in terms of a competition 
between two effects. At low temperatures the dominant 
dependence comes from the renormalization of the tun- 
nehng amplitudes Ti_2 in the presence of interaction Q. 
The critical current increases with temperature. Above 
the crossover temperature Tcross — (3/8)(2 — gp)vp/i:d 
the Josephson current decreases due to the shortening of 
the phase coherence length. Although the maximum is 
not very pronounced, the cross-over temperature shifts 
to higher values as the interaction strength increases 
(see Fig. This results in a wider temperature range 
in which the critical stays almost constant. It is evi- 
dent from Eq. (|2l| ) that the coefficient responsible for 
the anomalous temperature dependence vanishes in the 
absence of interaction, thus restoring the suppres- 
sion of the critical current (^0|) . For high temperatures 
T ^ Vp/2-Kd the suppression becomes exponential. 



4°)(T) cx T^/f" exp(-27rTd/t;F). 



(22) 



The full temperature dependence of the critical cur- 
rent, calculated by numerical integration of Eqs. ([T^), 
( |T8| ) is shown in Fig. ^. We see that for moderate strength 
of the interaction g ~ 1 the Josephson current will main- 
tain an appreciable value up to a temperature T ^ vp/d, 
which is of the order of 0.7 K for typical experimen- 
tal parameters vp = 3.0 10^ m/s and d = 3 /im [ p9[ . 
Moreover, the value of the critical current Ic°'\t = 0) « 
22 uA, (estimated for the parameters given above and 
Gi/{Ae'^) = 0.3) is large enough to be measured experi- 
mentally. 
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Note that we estimated the Josephson current assum- 
ing fixed Josephson phase difference between the su- 
perconductors. Thermal fluctuations of the Josephson 
phase would smear the critical current at temperatures 
T* ^ Ej = /i°^/2e, provided that the superconductors 
are coupled by the LL only. Using Eq. ( |l7|) for the non- 
interacting case (with Fi?\o) = 1), one obtains that 
the temperature T* is by a factor 27rGiG2/(4e^)^ ^ 1 
smaller than the characteristic temperature scale vp/d 
for the LL . Hence, in order to observe non-trivial tem- 
perature dependence of the critical current, one has to fix 
the phase difference between the superconductors, e.g., 
by means of an additional Josephson junction. 

b. The case vp/d^ A. 

In this limit, the electrons propagate fast and indepen- 
dently through the LL on a time scale 1/A. A typical 
contribution is depicted in Fig. ^b. The Cooperon (|lF 
can be approximated as 



Vl/L,+ (0,ri)*[,^(d,T2))(*L,_(0,T3)^'I_(d,T4)), 



n^"^(0,d;ri,...,T4) 

A 

where we substitute 

The constant C is determined by integration of the time- 
ordered single particle correlator of the LL, 

f3/2 

C= J rfT(T,*i,+ (0,r)*[^(d,0)). 

-0/2 

The temperature-dependent Josephson current is found 
to be 

4^)=/W(T)sinx, 
where the critical current is given by 

/W(T) = eA 

with 



_GlG2_ (6) 

(4e2/7r)2 ' 



1 


9p/4+l/3p-l 


'9p' 


5p/4+l/3p 


kpd 




. 2 . 





nd 
fivF 



££. + J_+l 



2 

TT 

n 



l3vF/2d 



dx sin 



Cp + C. 



cos\\{2t: d / j3vj) — cos2T:dx/ (3vf 



16 ^ 4gj 



The phase factor Q is given by 
Q = arctan 



Trdx \ , / Trd 
cot I - — tanh — - 
(Jvp J \pvj 



In the non- interacting case gp = = 2, F^'^ can be 
calculated explicitly: 



2 f 1 

- arctan . , , — r 

TT \smh(Tra/ pvp) 



At zero temperature we find i^i, = 1. The resulting 
Josephson current thus is independent of the distance d 
between the contacts, analogous to the result obtained in 
Ref . . At finite temperatures the Josephson current is 
suppressed. If T <C vp/d, the suppression is linear in T: 



i^c\t) 

#(0) 



ATd 

Vp 



In the interacting case, at T = 0, the critical current 
is suppressed, 



I'i'^O) cx I /{kpd) 



5p/4+l/Sp-l 



At finite temperatures, and for weak interactions we ob- 
tain 



i^'\t) 



3 

1 

2tt 



irTd 



Vp 



l/2+ffp/8+l/2gp 



2 — gp nTd 
(24) 



where we dropped terms of the order of 0[{2 — gp)'^ + {2 — 
gp){Td/vp)'^]. We find again an anomalous dependence 
on temperature, like the one we discussed in the case 
vp/d < A. 



V. DC-JOSEPHSON CURRENT THROUGH A 
RING 



In case of the ring with circumference L (Fig. |l|b), we 
should take into account the contribution to the Joseph- 
son current due to the topological part, see Eqs. (|4^), 
(^), and (H). The Cooperon for the ring in case of a 
symmetric setup d = L/2 ^ vp/A is evaluated along 
the same lines as before. It is then straightforward to get 
the Josephson current 



L (4e' 



e=±l 

X sin (x + enM„/2 + nJp/2)) jj^.j, 



(25) 



(23) where 
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kpL 



TrPvp/2L 



dx 



cosh(x) 

n 



■ cosh 



9p 



TrfjvF/2L 
2 



cosh(2x/gp) 



{Mp ~ 4/^)a; + ej^x 



X 

with 



l + E„exp(2^/3i;,nVi) 



1 + exp {2iT(3vjn'^ /L) cos{2nXj{x)) 



(26) 



Xj{x) — -arcos 



cosh ( 



The brackets, should be considered as the ther- 

mal average over the topological excitations weighted 
by the appropriate Boltzniann factor and subject to the 
topological constraints. 

For zero temperature this calculation involves only the 
ground state (see Ref. [|3) for details). The dependence 
of the critical current on /$ and has been found to 
show very rich behavior. In the present study, we will 
focus on the effect of finite temperatures. In particular, 
we will investigate how robust the structure, found in 
Ref. is against thermal fluctuations. 

Two remarks are in order at this point. First, through- 
out this section, we assume that the linearization point of 
the original electron spectrum corresponds to odd values 
of A^o (for even iVo the picture is the same, apart from a 
relative shift of /$ and /^). Second, the Josephson cur- 
rent depends periodically on both /$ and with period 
1. However, since the original problem has additional 
symmetries /$ — > — /$ and — (together with a 

change of sign of the corresponding topological numbers, 
Jj and Mj), it is enough to consider /$ and in the 
intervals < /$ < 1/2 and < < 1/2. 

It is instructive to discuss first the non-interacting case 
gp = 2, for r = 0. If fp + f<s> < 1/2 the ground-state con- 
figuration of Jj and Mj is found to be (Jp, Jg., Alp, M^) = 
(0,0,0,0); if fp + J^>l/2 the configuration is (2,0,2,0). 
Hence, the ground-state of the system can be changed by 
varying either the flux or the gate voltage. As a result, 
the Josephson current changes, see Eqs. (p5|), (p6|). This 
is illustrated in Fig. ^, where the critical current Ic (we 
write /,/ — Ic sin x where Ic can be positive or negative) 
is plotted as a function of and T at fixed /$ = 0.2. 
For T = 0, the critical current shows a maximum and a 
sharp jump at — 0.3 where the states (0,0,0,0) and 
(2,0,2,0) are degenerate. At this value of the gate volt- 
age, the number of electrons on the ring {Mp) increases 
by two. Since electronic states are doubly degenerate in 
spin and non-zero flux is applied, the two electrons will 
occupy the same (clockwise or counterclockwise moving) 
single-particle state. Therefore, the net current evpjp/ L 
increases by 2 quanta evp/L while the topological num- 
bers Mcr and Ja related to spin remain unchanged. At 



the jump, Ic changes sign. This reflects the fact that the 
ring acts as a 7r-junction {Ic < 0) in the state (2,0,2,0), as 
can be seen from Eq. (psf). Therefore, for non- interacting 
electrons, the critical current shows two jumps per period 
of the gate voltage dependence. The same is true for the 
dependence of Ic on the magnetic flux. 

This picture is correct for any generic point on the line 
U + U ^ 1/2. At the end points (/^,/$) = (0,0.5) 
and (0.5,0), no jumps of the critical current occur (one 
can say that two jumps in opposite directions merge to- 
gether). Instead, the critical current shows a resonance. 
The resonance occurs due to alignment of two spin degen- 
erate energy levels (for clockwise and counterclockwise 
moving electrons) with the chemical potential of the su- 
perconductors |2^ . 

At a flnite temperature, both the non-zero modes and 
the topological excitations are thermally activated. Ther- 
mal activation of the non-zero modes leads to an overall 
suppression of the critical current, as it has been dis- 
cussed for the wire in Section 4. Thermal activation of 
the states (0,0,0,0) and (2,0,2,0) will lead to a smearing 
of the jump. Moreover, at flnite temperature there will 
be a non- vanishing probability to activate other topolog- 
ical excitations which can contribute to the Josephson 
current. In the plotted temperature range, only one ad- 
ditional state (1,1,1,1) with one extra electron on the ring 
can be activated. As a result, the negative critical cur- 
rent of the state (2,0,2,0) (at fp > 0.3) will be partially 
compensated by the positive critical current of the state 
(1,1,1,1), the occupancy of which increases with temper- 
ature. Note that the jump at fp — 0.3 remains visible 
up to temperatures of the order oiT ^ vp/L ^ IK (for 
the parameters mentioned above and L — 2/im). Hence, 
the parity effect causing the jump is quite robust against 
thermal fluctuations. 

An important feature of the non-interacting case at 
T = is that the various possible ground-state conflgu- 
rations may differ by an even number of electrons only. 
The situation changes drastically when repulsive interac- 
tions are switched on. In addition to the states (0,0,0,0) 
and (2,0,2,0), the state (1,1,1,1) can act as a ground state 
conflguration This happens for parameters /$ 

within the range 1 + 3{gp/2f < 8[fp + (<?p/2)2/$] < 
3 -I- {gp/2)'^. Within this "strip", it is energetically 
more favorable to add a single electron, rather than a 
pair of electrons to the ring, due to repulsive electron- 
electron interactions. The Josephson current in the state 
(1,1,1,1) differs from the current in the states (0,0,0,0) 
and (2,0,2,0), see Eqs. (|25|), (|2|). For example, for 
gp = 1.75 and /$ — 0.2 (Fig. ||) the state (1,1,1,1) occurs 
in the range 0.259 < fp < 0.318. Indeed, one sees two 
pronounced jumps of Ic at the borders of this interval in 
Fig. ^ (at low temperatures). Generally, for interacting 
electrons the critical current shows four jumps per period 
of the gate voltage dependence . 

Similar jumps are seen also at the dependence of the 
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critical current on the flux, Fig. |[ Moreover, Ic is a 
stepwise function of for T = (this can be seen from 
Eqs. (|2^), (|6|); the flux /$ enters to these equations only 
implicitly, via topological numbers). Depending on the 
gate voltage, the critical current can show zero, two, or 
four jumps per period of the flux dependence | |2^ . 

The state (1,1,1,1) is the ground state in a strip of 
width Sff, = [1 - (5p/2)2]/4. This determines the energy 
SE ~ [1 — {gp/2)'^]TTVp/{gpL) needed to create topologi- 
cal excitations. The features related to the state (1,1,1,1) 
will be smeared at temperatures T ~ SE. Therefore, for 
weak interaction 2 — gp ^ 1, the interaction effects will 
dissapear at much lower temperatures T <C vp/L than 
the parity effects. For example, the features related to 
the conflguration (1,1,1,1) in Figs. ||, || are seen only in 
the temperature range T < 6E ~ 0.1 K whereas the over- 
all dependence is robust up to the temperatures T ~ IK. 
However, it is worthwhile to stress that the (1,1,1,1) state 
survives at much higher temperatures when the interac- 
tion strength is increased. 

The behavior we described here is rather generic for 
all values of fp, /$ and for various values of the inter- 
action. What is specific is the configuration of the two 
superconductors: they are connected symmetrically to 
the ring. If the points on the ring at which the elec- 
trodes are attached would form a generic angle, a more 
complicated interference pattern would arise. In the sym- 
metic setup, the maximum Josephson current occurs ei- 
ther at X = or at X = In the nonsymmetric setup 
the maximum Josephson current would occur at a value 
x( Jp, J<T, Mp, Ma) which depends on the values of topo- 
logical numbers. The critical current should then be 
found by maximizing the resulting function of the phase 
difference. 



average of the corresponding Heisenberg operator. Us- 
ing the interaction representation with the unperturbed 
Hamiltonian Hq, see Eq. (^), one obtains 

m = {u\t)i{t)u{t)), 



U{t)=Tcxp{-i HT{t')dt'). 



(27) 



We expand ( p7| ) to the fourth order in Ht and keep the 
Josephson terms in the current. These are proportional 
to exp(±2jeyt). As a result the Josephson current is 
given by an expression which has the same structure as 
Eq. ( p^ . The integrals are now taken over real times. It 
is convenient to depict the times t,ii,i2,^3 of tunneling 
events on the Keldysh contour |Q (the Josephson cur- 
rent is calculated at a time t). Again, we will consider 
two cases of long {vp/d <C A) and short {vp/d » A) dis- 
tance between the contacts. The relevant diagrams are 
shown in Fig. ^ for both cases. We restrict our consider- 
ation to the case of a quantum wire at zero temperature. 

a. The case vp/d <C A. 

For a large distance between the contacts the tunnel- 
ing of two electrons to/from a superconductor is a fast 
process on the time scale of their propagation through 
LL. The Josephson current is described by diagrams of 
the type shown in Fig. 0a. The Josephson current is then 
given by 

2 2 ^102 



Ij{t)=ATr^ev'p^^ 



x3?e 



±2ieVt 



(4e2)^ 
dt'e=F*«^*'n(t') 



(28) 



VI. AC-JOSEPHSON EFFECT 

The effect of a finite DC bias voltage eV ^ 2A applied 
between the superconductors SI and S2, will be twofold. 
First of all, the phase difference x between SI and S2 will 
acquire a time-dependence, according to the Josephson 
relation x = wj = 2eV. As a result, the Josephson 
current will oscillate as a function of time at a frequency 
ujj (AC- Josephson effect, see Ref. [|o|). Secondly, a DC 
subgap current will be induced, due to Andreev reflection 
at both junctions. This current is dissipative, energy will 
be dissipated in the LL. In a typical experiment one thus 
will find a current with both a DC and an AC component. 
In this section, we will mainly concentrate on the AC 
Josephson current, and estimate the DC component at 
the end. 

In the presence of a bias voltage V between the super- 
conductors, the imaginary time formalism cannot be ap- 
plied and Josephson current is found by calculating the 



where Tl{t) = 11^'*^ (0, c?; it, it, 0, 0) is the Cooperon prop- 
agator (^^ in real time taken at coinciding time argu- 
ments. The leading contribution stems from the terms in 
Eq. d) with n = ±1, 

n(t) = 

2po Y[{[l + ikp{vjt + d)][l + ikp{vjt-d)]}-^/s,^ (29) 

In particular, for non-interacting electrons {gp — 2) we 
obtain 



27revp G1G2 . eVd 
— ; — . . sm I 2eKt 



d (4e2)2 



Vp 



(30) 



This result means that the Josephson current acquires 
an additional phase shift due to the propagation of elec- 
trons between the contacts. For interacting electrons we 
computed the Josephson current numerically. We split 
I J into sinusoidal and cosinusoidal components, 
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T (f) - 27ret;_F G1G2 / 

d (4e2)2 \kFd 



2/Sp-l 



X <! Js (-^,gp ] sm{2eVt) 



+ Jc 



eV 



vp/d 



, gp cos {2eVt) 



(31) 



The amplitudes Jc{s) of the two components and the 
phase if = — arctan( Jc/ Js) of the Josephson current 
are shown in Fig. ^ as functions of the voltage for two 
values of the interaction parameter, gp = 1.75 and 1. 
One sees that the deviation from the simple result ( |30| ) 
for non-interacting electrons (which corresponds to Js — 
cos{eVd/vp), Jc — — sin^eVd/vp), and ip — eVd/vp) 
increases with the increase of the interaction. This de- 
viation becomes striking in the dependence of the AC 
current amplitude Ja = \/ Jg + Jc on the voltage, Fig. ||. 
Apart from the non-interacting case {Ja(y) = const), 
one sees pronounced oscillations of the current ampli- 
tude. These oscillations are due to the difference in the 
velocities of the charge (vp) and spin (va-) excitations. 
The period SV of the oscillations corresponds to 27r dif- 
ference between the phases of charge {eVd/vp) and spin 
{eVd/vp) excitations. Using the relation Vp = Ivpjgp 
we obtain e5V/{vp/d) = 2ti{\ - gp/2)-^ ~ 50.4, 12.6 for 
Qp = 1.75 and 1, respectively. This is in very good argee- 
ment with the period of oscillations in Fig. ||. Therefore, 
the AC Josephson effect can be used as a tool for the 
observation of spin- charge separation in the LL. 

b. The case vp/d^ A. 

At short distances between the contacts, the two elec- 
trons propagate fast through the LL on the time scale 
1/A. The relevant diagrams are similar to the graph 
shown in Fig. |^b. The main contribution to the Joseph- 
son current comes from the integration of the two-particle 
propagators of the type n^^'(0, d; it, iti,it2, its), Eq. ( |l5| ) 
(with possible permutations of creation and annihilation 
operators) over the range t — ti ^ t2 — t^ ^ d/vp and 
t—t2 ^ 1/A (see Fig. ^p). For this reason, we can present 
n^^' as a product of two single-particle propagators and 
integrate the latter over the "fast" variables t — ti and 
t2 — ts (from to 00) as we did for the DC case. The last 
integration over the "slow" variable s — t—t2 involves the 
product of two anomalous Green functions with an ex- 
ponent containing the time-dependent Josephson phase: 

/•OO 

/ dsF+^{Q,is)Fs2{0,is)exp{ieVs). 
Jo 

Hence, for short distances between the contacts the pres- 
ence of LL does not influence the voltage dependence of 
AC Josephson current. The latter is still given by the 
simple formula 

Ij{t) = (2/7r)X(ey/2A)4^)(0) sin(2eyt) , 



where K{x) is an elliptic integral and Ic''\o) is the crit- 
ical current in the DC case (cf. Eq. ( p3|) in the limit of 
zero temperature) . The effect of the interaction is only to 
reduce the value of the critical current while its voltage 
dependence is analogous to that of the critical current of 
a conventional Josephson junction |^^ . 

We conclude this section with an estimate of the dis- 
sipative DC current due to Andreev reflection at both 
junctions. For a single junction between a superconduc- 
tor and a LL with repulsive interactions, the subgap cur- 
rent Is{V) as a function of the apphed voltage V is given 
by ||,||,|| Is(y) - y|y|2/9p-i. For the system of 
Fig. la, which consists of two junctions in series, the 
lowest order contribution to the subgap current stems 
from sequential tunneling. Employing a rate equation 
approach, we find for this contribution 



IsiV) = 27re-SiS:2ey 



■(4e2)2 



2eV 



vpk[ 



(2/5, 



2(GiG2) 



r(i 

Sp/2 



2/<?p) 

2/ffp 



(G2)9p/2 + (G2)9p/2 



(32) 



Comparing this result with the critical current we see 
that the dissipative component is much smaller at low 
voltages. In order to get a complete description at finite 
voltages, one has to solve the corresponding equation for 
non-linear RSJ-modcl pl. This will be discussed in a 
forthcoming publication |34| . 



VII. DISCUSSION 

In this paper, we studied the AC and DC Josephson 
effect in a single mode quantum wire and quantum ring 
connected to two superconductors by tunnel junctions. 
Repulsive interactions were treated in the framework of 
the Luttinger model. Interactions were found to have a 
drastic influence on both DC and AC Josephson effect. 

The critical current is suppressed by interactions at 
zero temperature. The results depend on the ratio be- 
tween the characteristic energy fivp/d of the ID electron 
system and the superconducting energy gap A. For large 
distances between the contacts d ^ hvp/A in the pres- 
ence of interactions, there is a competition between ther- 
mal suppression of coherent two-particle propagation in 
the wire and activation of tunneling at the junctions at 
low temperatures. As a result, the critical current shows 
maximum as a function of temperature. At even higher 
temperatures, fc^T ^ hvp/2'nd, the suppression becomes 
exponential. 

In our model it is assumed that the superconducting 
electrodes do not influence the uniformity of the potential 
along the quantum wire, since they are separated from 
the wire by thick barriers. It was argued in Ref. 13^ that 
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a non-uniform potential in the wire will lead to an effec- 
tive change of the boundary conditions for the electronic 
wave function, which in turn could strongly affect our 



results. However, a recent calculation |24| of the Joseph- 
son current through an interacting quantum wire of finite 
length is in agreement with our results. This indicates 
that the results obtained are robust with respect to the 
specific choice of boundary conditions. They rather de- 
scribe generic properties of the superconductor-Luttinger 
liquid system. 

If a finite voltage V is applied between the junctions, 
the AC Josephson effect occurs. The AC current ac- 
quires phase shift proportional to the distance between 
the tunnel junctions. Moreover, the amplitude of AC cur- 
rent depends on the voltage in an oscillatory fashion due 
to spin-charge separation. The corresponding period de- 
pends on the ratio of the velocities of the spin and charge 
excitations in the LL. 

A quantum wire closed to a loop (or quantum ring) 
shows interesting parity effects. Boundary conditions on 
the electronic wave functions result in a discrete set of 
quantum numbers, related to the number of particles and 
angular momentum. We showed how these numbers can 
be tuned by applying a gate voltage and a magnetic flux, 
and calculated the corresponding dependence of the criti- 
cal current on these parameters. This dependence shows 
a rich behavior, which can be detected in an interfer- 
ence experiment employing a SQUID. We showed that 
the dependence is robust to thermal fluctuations up to 
experimentally measurable temperatures. 
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FIG. 6. The critical current through the ring at a fixed 
value of the gate voltage = 0.2 as a function of the 

flux and the temperature (i = TL/nvp) in the interact- 
ing case Qp = 1.75 (the critical current is normalized to 
2nevF G1G2 s 
~^^{4eyhf'' 

FIG. 7. Relevant diagrams for AC Josephson current: (a) 
vp/d^ A and (b) vp/d <^ A. The shaded area indicates the 
electron-electron interaction. 

FIG. 8. The voltage dependence of sinusoidal (solid line) 
and cosinusoidal (dashed line) components and the phase 
(dotted line) of AC Josephson current; (a) gp = 1.75 and 
(b) gp = 1. 

FIG. 9. The voltage dependence of the amplitude of AC 
Josephson current. Here, gp — 2, 1.75, 1, 0.75, 0.5, 0.25 for the 
curves from top to bottom at zero voltage. 



FIG. 1. The geometries discussed in the text: (a) 

one-dimensional wire connected to two superconductors by 
tunnel junctions. The distance between the junctions is d. (b) 
Ring with circumference L connected to two superconductors 
by tunnel junctions. The distance between the junctions is 
L/2, the ring is threaded by a magnetic flux #. 



FIG. 2. Relevant diagram for Josephson tunneling in the 
limiting cases (a) vp/d A and (b) vp/d 3> A. The shaded 
area indicates the electron-electron interaction. 



FIG. 3. The critical current of the wire as a function of the 
temperature {t = Td/2nvF) for various values of the interac- 
tion strength gp = 1.0, 1.25, 1.75, 2.0. 



FIG. 4. The critical current through the ring (normalized 

to ^^^^^ , Pt'^^.o ) at a flxed value of the flux = 0.2 

L {Ae-^/ny 

as a function of the gate voltage and the temperature 
{t = TL/nvp) in the non-interacting case. 



FIG. 5. The same as in the previous figure for the interact- 
ing case; for this plot we choose gp — 1.75 but the result is 
rather generic for repulsive interaction (the critical current is 

,. J ^ 27re«F G1G2 , 
normalized to — - — „ ,-^s, )- 
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